Role of rotational symmetry in the magnetism of a multiorbital model 
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Effect of rotationally-invariant Hund's rule coupling on a magnetism of multiorbital Hubbard 
models is studied within a dynamical mean field theory framework. Comparison of static magnetic 
susceptibilities and local densities of states of two- and three-orbital models of a complete rotationally 
invariant Coulomb interaction and a "density-density" Hartree type interaction shows the different 
role of spin-flip interactions for different band fillings. In the particle-hole symmetric case the 
Mott-Hubbard physics dominates due to the strong effective Coulomb interaction, while for the 
multiple electronic configurations away from half-filling (two electrons in the three band model) 
the formation of local magnetic moments due to Hund's exchange interaction becomes the most 
significant effect for itinerant magnetic systems. A shift of the temperature of magnetic ordering 
due to the rotationally-invariant Hund's rule coupling is found to be the largest in a three-orbital 
model with a two-electron occupancy where the single particle spectrum is metallic and is not 
sensitive to different forms of the Coulomb vertex. A larger enhancement of the effective mass in a 
model with a rotationally-invariant interaction is discussed. In the half-filled case we find a drastic 
change in the density of states close to the Mott transition which is related to the spin-flip Kondo 
fluctuations in a degenerate orbital case, while the corresponding shift of the magnetic transition 
temperature is relatively small. It is shown that a change in the ground state degeneracy due to a 
different symmetry of the Coulomb interaction in the density-density model leads to a breakdown 
of the quasiparticle peak at the Fermi level in the proximity of a Mott transition on the metallic 
side. We discuss the relevance of rotationally-invariant Hund's interaction in the transition metal 
magnetism. 



PACS numbers: 71.27.+a, 75.10.Jm, 71.10.Fd 
I. INTRODUCTION 

In many materials with strong electronic correlations 
spin and orbital fluctuations play an important rol&U Ex- 
perimental spectroscopic datcP and analysis of a simple 
modeP show that the local Coulomb interaction of va- 
lence electrons is sufficiently screened, while the higher- 
order Slater integrals are almost unaffected by a metallic 
screening. This gives rise to important effects, associ- 
ated with the Hund's exchange coupling J in strongly- 
correlated magnetic compounds. In the case of relatively 
large onsite Coulomb repulsion in transition metal com- 
pounds a ferromagnetic state can be stabilized by orbital 
ordering of the orbitals in the valence shelPP A vari- 
ation of an effective hybridization between different or- 
bitals can lead to an orbital-selective Mott transitiorP^ 
which was observed experimentally in La„+iNi n 03 n +r^I. 
In the strongly correlated metals away from the particle- 
hole symmetry the effective Coulomb interaction is well 
screened and the relevant electron interaction is related 
to a Hund's rule coupling with complicated multiplct 
fluctuations-^. The physics of such Hund's metals is dom- 
inated by the spin-flip type of fluctuations which lead to 
a metallic state with correlated local-moments and a re- 



duced spectral weighlP^. 

The effective way of making a non-perturbative de- 
scription of strongly correlated electronic systems is a 
self-consistent mapping of a complex lattice problem to a 
single impurity with an energy-dependent external bath. 
In this dynamical mean field theory (DMFTp^l the re- 
sulting Anderson impurity problem is solved numerically. 
Until recently numerical solvers for this problem could 
treat the Hartree-type density-density Coulomb interac- 
tion, while an inclusion of the off-diagonal terms l ed to 
an exponential increase of the computational tim d 14 * 15 -*. 
A significant progress has been achieved with the ap- 
pearance of the continuous-time quantum Monte-Carlo 
algorithms^. An application of the algorithm to a single 
Anderson impurity problem with a 3d valence shell (Co 
atom on a metallic surface) and taking into account a full 
multiorbital character of the Coulomb interaction results 
in a shift of the Kondo peak in the density of states and 
a reduction of the Kondo temperature^*^. 

A realistic description of the strongly-correlated elec- 
tron materials such as transition metal systems with par- 
tially filled 3d bands is often impossible without taking 
into account multiplet effects in the valence shell. In- 
deed previous LDA+DMFT calculations for the density- 
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density type interaction overestimate the Curie temper- 
ature of irorP^ and various manganites^ by a factor of 
1.5 — 2. This inconsistency can be associated either with 
the absence of non-local correlations in the dynamical 
mean- field scheme^ or with the lack of spin- flip and pair- 
hopping terms which are prese nt in the realistic multior- 
bital Hubbard-like modeP2EH I n order to separate these 
contributions we study the Hubbard model in the limit of 
infinite dimensions w here t he local dynamical mean field 
theory becomes exadP^S. 

In recent years an extensive study of the two- and 
three-orbital Hubbard models in infinite di mensio ns 
with Numerical Renormalization group 
self-energy functional approaches equation of mo- 
tion method^, highly optimized Hirsch-FydJ^, exact 
diagonalizatiorPS and continuous time Quantum Monte 
Carlo (CT-QMC^pHaU so lvers has been performed and 
the phase diagrams of the models are established. In- 
sulating phases occur at integer electronic fillings. In 
the half-filled case the antiferromagnetic behavior is ob- 
served for any U larger than the critical value of the 
metal-insulator phase transition, which is of the order of 
bandwidth. Away from particle-symmetric regime of the 
model a ferromagnetism coupled with the orbital order 
of the electrons in the lattice can be stabilized at large 
values of U, otherwise a paramagnet or in the case of 
filling with 2 electrons in a 3-orbital shell the frozen lo- 
cal momentPS phase is observed. An application of the 
model to various transition metal oxides with a large 
Hund's coupling reveal a strong enhancement of the ef- 
fective mass depending on the valence band electronic 
fillingfP. 

Despite of a general understanding of the low- 
temperature phase diagram of the multiorbital Hubbard 
model, a careful comparison with the simplified density- 
density Coulomb interaction model, which is generally 
used in applications to real materials, needs to be per- 
formed. One of the experimentally accessible observable 
quantities is the critical temperature of magnetic phase 
transitions. In our recent letter [37] the static magnetic 
susceptibility of the half-filled two-orbital paramagnetic 
Hubbard model was obtained using the continuous-time 
interaction-expansion quantum Monte-Carlo solver^. 

In this paper we discuss the dependence of the mag- 
netic phase transition on the electronic filling in two- and 
three-orbital Hubbard models for the general rotation- 
ally invariant and density-density Coulomb interactions. 
The interaction expansion v ariant of the continuous-time 
quantum Monte-Carlo solver^^U is used for the calcula- 
tions. The paper is organized as follows. In section II 
the general multiorbital rotationally invariant Hubbard- 
like model is introduced. In section III the results for the 
static magnetic susceptibilities for two- and three- orbital 
models on an infinite-dimensional Bcthc lattice for differ- 
ent electronic fillings are presented. We discuss the im- 
portance of spin-flip Hund's rule effects in the half-filled 
model in section IV and away from half-filling in section 
V. In appendix [X] we show that the reduced degeneracy 



of the ground state of the "density-density" model leads 
to the disappearance of a quasiparticle peak at the Fermi 
level. This behavior is effectively traced in a hybridiza- 
tion expansion in auxiliary dual fermions, which uses ver- 
tex functions of the "atomi c" pro blem as a building block 
for the perturbation theorjE^HSB. 



II. MODEL 

The multiorbital extension of the Hubbard model can 
be written in the following form: 



(ij)acr 
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where t denotes the (orbital independent) hopping be- 
tween adjacent sites i and j of the lattice and sets the 
energetic scale of the model, a and a are orbital and spin 
indices, is the chemical potential and H lnt is the local 
onsite Coulomb interaction. Without symmetry break- 
ing the Hamiltonian of the onsite Coulomb interaction 
obeys the rotational invariance and the conservation of 
spin and can be expressed in terms of the squared total 
spin S 2 , squared angular momentum L 2 and number of 
particles N operators. For the shell with conserved I = 1 
angular momentum and s = 1/2 spin of each electron it 
reads^ 
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Being expressed in the basis of Fock states the Hamilto- 
nian ([2| takes the form introduced in Ref. [15] and can be 
split into a part which corresponds to a density-density 
type of the interaction and the remaining off-diagonal 
terms, which include spin-flip and pair-hopping terms: 
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Here a denotes the opposite of a. The model |3| keeps its 
degeneracy in a cubically symmetric environment and as 
such (with an introduction of parameter U' , here equal 
to U — 2 J) is used for a description of a three-orbital t 2g 
and also two-orbital e g electronic shelld^, obtained for 
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Table I. Eigenfunctions of the rotationally invariant (top) and 
"density-density" (bottom) Hamiltonians, which belong to 
the subspace of the double electronic occupancy for the two- 
orbital model, n denotes number of electrons, E — energy 
of the eigenstate, S — total spin, Sz - spin projection on the 
z-axis, g — degree of degeneracy. The ground state is typeset 
in bold. 
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Table II. Eigenfunctions of the rotationally invariant (top) 
and "density-density" (bottom) Hamiltonians, which belong 
to the subspaces of two- and three- electrons occupancy for 
the three-orbital model Q. The notations are the same as in 
Table [I] apart that only the eigenstates nonequivalent under 
spin and orbital space rotations are included in the table (the 
number of the other equivalent states is denoted with "st"). 
The ground states for each electronic occupancy are typeset 
in bold. 



crystal field splitting of the d-shell. Note that equation 
([2]) is not applicable for the two-orbital case. 

Several remarks should be done about the single- 
atomic limit of the model. At t = rotationally invariant 
and "density-density" models have different degeneracy 
of the ground states in various electronic occupation sec- 
tors. The spectrum of the model is shown in tables [I] 
and[IT]for two- and three- orbital complexity respectively. 
The tables are analogous to the ones in Ref. [3H and sup- 
plementary material of Ref. 11, These are exactly the 
pair-hopping and spin-flip terms which entangle the Fock 
states in ^ and make the spin triplet a ground state in a 
half-filled rotationally invariant model. As it was pointed 
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Figure 1. (Color online) The temperature dependence of the 
inverse magnetic susceptibility of a single unhybridized atomic 
orbital, defined by full rotationally invariant Hamiltonian Q 
(solid lines) and Q (dashed lines) in the particle-hole sym- 
metric regime at U = 3.6, J = 1 for orbital dimensionalities 
1,2,3. 



out by Pruschkc and BullctM' the change in the character 
of the ground states has a profound impact on the exci- 
tation spectra. Kondo-type spin-flip processes are clearly 
possible in the atom with a multiplet classification of mul- 
tielectron levels, where S z can be changed to S z ± 1, pre- 
serving the atomic energy. Hund's rule guarantees that 
the ground state in a rotationally invariant model has the 
largest possible (S* 2 ), so that the processes S z — > S z ± 1 
are indeed allowed and the contribution to the central 
quasiparticle peak of Kondo scattering is present. 

In the "density-density" case the ground state has a 
lower degree of degeneracy and consists only of Fock 
states with all spins aligned in one direction which cor- 
responds to the Ising-type interaction. The spin projec- 
tions of different ground states differ by AS Z = 1 for the 
2-band model and by AS Z = 3/2 for the 3-band model, 
so no scattering process in the impurity can induce a gap- 
less Kondo excitation. This behavior is observed in the 
excitation spectra of the lattice models in the vicinity of 
a Mott transition and is discussed in section IV and its 
analytic treatment in this case is presented in Appendix 

El 

The change in the degree of degeneracy of the ground 
state of the "atomic" solution of the model has a clear 
consequence on the Curie-law dependence of inverse mag- 
netic susceptibility on temperature, shown in Fig. [T] In 
the low-temperature regime T <C J < U the Curie con- 
stant is obtained as a response of a single S = 1 or 3/2 
spin for two- and three-orbital models respectively. The 
density-density type Hamiltonian possesses only states 
with \S Z \ — S and thus has larger Curie constant and 
magnetic susceptibility. Indeed at T < J, = 2/3, 

Cat= 2 = !> C N=3 = 5 / 4 > C N n =3 = 9 / 5 - In thc high- 
temperature regime weights of all states become the same 

and the Curie constant is defined solely by the dimension- 
ality of the Hilbert space and has no difference between 
density-density and rotationally invariant models. 
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III. STATIC MAGNETIC SUSCEPTIBILITY 
CALCULATION 

We proceed now with the lattice calculations. In the 
limit of infinite dimensions the linked-cluster interaction 
expansion^ of the Hubbard model is reduced to the anal- 
ogous of the Anderson impurity modeP^, with the effec- 
tive hybridization function A arT (iui n ) equal to 

A acr (iuj n ) = t 2 g ag -(iuj n ), (5) 

where the spin-polarization can be imposed with a index 
to stabilize the anti-ferromagnetic ordering. The action 
for the lattice problem then can be replaced with the 
following: 

gimp _ _ y^y^ (iu n A aCT (iw n ))c* nQ , CT c w „ ao .+ 

+ J dTH int [c,c*}-2 J dThS z (r), (6) 
o 

where oj n are Matsubara frequencies, S z = \ X) a (?W — 
n a \) is the total spin projection along the z axis and 
electronic magnetic moment is taken as unity. In order to 
obtain the static magnetic susceptibility a small external 
magnetic field h is applied along z axis. The magnetic 
susceptibility then reads: 




Figure 2. (Color online) Curie-Weiss dependence of the in- 
verse magnetic susceptibility x ^^ 1 ) over temperature of the 
half-filled two- and three-orbital Hubbard model at U/t = 
4, J = 1.2t and U = 8t, J = 1.2t (for the two-orbital case 
only). The curves for the rotationally invariant interaction 
are named "full" and plotted with straight lines and filled 
points, while the curves which correspond to the density- 
density type of interaction are named "nn" and plotted with 
dashed lines and empty points. The dotted lines are linear 
fits of the inverse susceptibilities in the paramagnetic regime. 
The Neel temperature is defined as the crossing point of in- 
terpolation lines with the X -1 ^) — axis. The temperatures 
obtained are T^=2 = 0-42 in the rotationally invariant model 
and Tf^L 2 = 0.49 in the density-density variant of the two- 
orbital model. In the three-orbital case the corresponding 
values are I$£ 3 = 0.7 and T^L 3 = 0.83. 
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where q = 0, X are generalized Brilloin zone points 
corresponding to ferromagnetic (0) or antiferromagnetic 
(X) staggered magnetic field. Finally, the critical 
(Curie/Neel) temperature is obtained by a linear inter- 
polation of the Curie- Weiss law: 



C 



X = 



T - T, ' 



(8) 



where C is the Curie constant. Since at the limit of in- 
finite dimensions the lattice effects are included at the 
dynamical mean-field level, the Weiss temperature ob- 
tained in the aforementioned way is indeed the magnetic 
ordering temperature^. Note that it is also positive for 
the antiferromagnetic ordering as a response of bipartite 
(pil) lattice on the staggered field ([7| applied with q = X. 
This corresponds to the x(<f = X) susceptibility in finite- 
dimensional lattices. 

The results for the inverse of the static magnetic sus- 
ceptibility as a function of temperature for the models 
with full rotational invariance and density-density type 
of the interaction in the particle-hole symmetric case are 
plotted on Fig. [2j First, it is seen that exactly as in the 
case of an "atomic solution" on Fig. [T]an increase of the 
total number of orbitals leads to an enhancement of the 
static magnetic susceptibility in the half-filled case, yet in 



the lattice it is also followed by a corresponding increase 
of the Neel temperature. The ratio between T c in half- 
filled rotationally-invariant two and three orbital models 
is equal to 1.7 which is a bit less than the 15/8 ratio of the 
mean field critical temperatures in the Heisenberg model, 
for which the magnetic ordering temperature in the mean 
field approximation is defined as T c = J'zS(S+l)/3 with 
the magnetic exchange constant J and z nearest neigh- 
bors in the lattice. The difference between these values 
is associated with the charge fluctuations absent in the 
Heisenberg treatment. Comparing curves for the two- 
orbital model at U/t — 8 and U/t = 4 with the same 
value of J = 1.2t reveals that an increase of Hubbard U 
and thus a reduction of the effect of charge fluctuations 
leads to a small change in the Curie constant and does 
not affect the Neel temperature in this regime. 

The important consequence of taking into account 
spin-flip and pair-hopping terms is that it leads to a de- 
crease of a magnetic ordering temperature. It is clear 
that the spin-flip terms generally help to rotate large lo- 
cal magnetic moments and reduce the critical temper- 
ature. This behavior has a prototype in a comparison 
of critical temperatures of Ising and Heisenberg models. 
In the half-filled case on Fig. [2] the magnitude of this ef- 
fect is approximately 20% both for two- and three-orbital 
models. 

The curves for the inverse of the static magnetic sus- 
ceptibility as a function of temperature away from half- 
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Figure 3. (Color online) Curie- Weiss dependence of the in- 
verse magnetic susceptibility x~ l {T) over temperature of the 
two-orbital model at the single-electron filling and the three- 
orbital model at two-electron filling. The curves are obtained 
at U — 8, J — 1.2t. The notations are the same as in Fig. 
[2] The obtained magnetic transition temperatures in the 
three-orbital model at one-third filling differ by a factor of 
2 and have values T^5 3 . n _ 2 = 0.14 in the rotationally in- 
variant model and T^L 3 . n=2 = 0.27 in the "density-density" 
model. 



filling arc plotted on Fig. [3j The case of the one-third 
filling (with n = 2 electrons) in the three orbital model is 
of particular interest. Previously it has been shown that 
a state with a frozen moment is realized at low temper- 
atures at moderate values of t7^, whereas an increase of 
U leads to the ferromagnetic spin ordering^. We thus 
impose a self-consistency condition ^ without spin in- 
version in this case^l. It is first seen, that the static 
uniform magnetic susceptibility is larger than in the half- 
filled case. An explanation for this fact is that the degree 
of degeneracy of the ground state and thus the magnetic 
moment of the shell with two electrons is larger than of 
the half-filled one by a combinatorial factor (see Table 
|fl| ). Larger magnetic susceptibility in the one-third filled 
three-orbital shell leads to an enhanced ratio between 
Curie temperatures of the density-density and rotation- 
ally invariant Coulomb interaction models which in the 
described case reaches a factor of two. This enhancement 
is related to the fact that correlation effects are almost 
completely dominated by the Hund's parameter J, while 
the effective Coulomb interaction is sufficiently reduced. 
The difference between the half-filled and one-third filled 
cases is covered in greater detail in Discussion. 

A large doping of the system which corresponds to 
n = 1 average electrons number per site in the two or- 
bital model leads to the obvious result. It is seen on 
Fig. [3] that in this case no spin-flip or pair-hopping pro- 
cesses which involve two electrons can occur, and thus 
there is no difference between rotationally invariant and 
density-density models. The fact that the line of the in- 
verse magnetic susceptibility crosses at a negative value 
corresponds to the fact that the model in the described 
regime stays paramagnetic at zero temperature^. 
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Figure 4. (a) Family of densities of states of the two-orbital 
Hubbard model for various values of U in the vicinity of 
the metal-insulator transition, J = 17/4, T = 0.04t. All 
curves are produced by using the analytic continuation by 
the maximum entropy methocP. The upper part of the 
graph contains results for the rotationally-invariant Coulomb 
interaction the lower — for its density-density part Q. 
(b) Density of states at the Fermi level A(0) for rotation- 
ally invariant and density-density models as a function of U, 
J = U/4 at temperatures 0.04t, 0.05t, 0.07t. The points were 
obtained from the extrapolation of Matsubara Green's func- 
tion to the zero value and showed no difference from the value 
PG(t = P/2)=A(u = Q)\ t ^o. 



IV. DENSITY OF STATES 

In this section we study the low-temperature behav- 
ior of rotationally invariant and density-density models 
and probe the Mott transition by changing the onsite in- 
teraction parameter U in two and three-orbital models 
at various fixed temperatures at half-filling. In order for 
the Mott transition not to get hidden by the antiferro- 
magnetism the paramagnetic model is studied and the 
symmetry against spin inversion in the effective impurity 
model is enforced. The external magnetic field is kept at 
zero value. 

Densities of states for the two-orbital model in the 
vicinity of a Mott transition obtained by analytic 
continuation^ for various values of U are plotted on Fig. 
[4^,. The parameters are chosen to be the same as in 
Fig. 7 of Ref. E] and J is fixed at (7/4, yet the finite- 
temperature behavior at T = 0.04t is described here. 
The curves for the rotationally invariant interaction |3]) 
are reminiscent of the single-band Hubbard model and at 
U < U c ~ 3.3i show a quasiparticle peak at the Fermi- 
level, coupled with shoulders from the Coulomb interac- 
tion at ±{7/2, whereas at U > U c the charge gap in the 
density of states is opened. The finite temperature effects 



Figure 5. Family of densities of states of the two-orbital Hub- 
bard model for various values of temperature for the mod- 
els with the (a) rotationally-invariant Coulomb interaction 
at U = 3.14t, J — (7/4, (b) density-density interaction at 
U = 2M, J = U/4. 



lead to the lack of "pinning" of the central quasiparticle 
peak and thus its height is not fixed for different values of 
U in the metallic side. This effect is clearly seen on Fig. 
|4Jd, where the value of density of states at Fermi level 
is plotted as a function of U at different temperatures. 
Here a decrease of the temperature leads to "straighten- 
ing" of the line of A(u> = 0) values in the metallic part 
and so one expects the "uncorrected" value of DOS to 
be restored at zero temperature. 

A comparison of the rotationally invariant and 
"density-density" curves at Fig. and |4]d reveals the 
same effect previously observed at zero temperature^, 
namely that taking into account only "density-density" 
part of the interaction leads to a breakdown of the central 
quasiparticle peak at Fermi level and thus to a sufficient 
reduction of the critical U to « 2.9t. In order to explore 
the difference between models in greater detail we fix the 
value of U and plot several densities of states at U = 3.14i 
and U = 2.9t for rotationally-invariant and "density- 
density" models respectively on Fig. [5] These points are 
chosen close to U c so that at the temperature T = O.lt 
both models show no quasiparticle peak at Fermi level. 
With a decrease of temperature to 0.035i the three-peak 
structure is observed for the rotationally-invariant model 
while the gap persists in the "density-density" case. The 
fact that the central peak emerges with decreasing of the 
temperature reveals its truly Kondo nature. This also in- 
dicates that absence of this peak is related to the lack of 
Kondo-type resonances in the "density-density" model. 

We conclude this part by comparing densities of states 
of the half-filled three-orbital model for the rotationally 
invariant and "density-density" interactions at various U 
in the vicinity of a Mott transition on Fig. [| The curves 
are almost identical to the two-orbital ones and the same 
behavior of the central peak is seen. The main difference 
is that the value of critical U is even more suppressed 
to w 2.6t and rs 2.2t in the rotationally invariant and 
"density-density" models respectively. This reduction of 
critical U with the number of orbitals directly follows 



Figure 6. Family of the density of states of the three-orbital 
Hubbard model at half-filling for various values of U in the 
vicinity of the metal-insulator transition, J = U/4, T — 0.05t. 
The curves are plotted in the same fashion as on Fig. |4j upper 
part corresponds to rotationally invariant interaction, lower - 
to its "density-density" part. 
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Figure 7. Densities of states of the one-third filled (2 
electrons) three-orbital model at U = 2.6t, J = (7/4 for 
T — 0.04 and 0.05t obtained by maximum entropy analytic 
continuation^. 

from the formula (|3| . Keeping J = restores the U c = 4i 
value of the critical interaction strength^. 

A more detailed discussion on the mechanism of forma- 
tion of the central quasiparticle peak is given in Appendix 

m 



V. DISCUSSION 

The results obtained demonstrate a drastic difference 
between the two different regimes, namely for the half- 
filled and the one-third filled case. The physics of a 
particle-hole symmetric system is clearly of the Mott 
type. Hubbard bands are formed near the point of metal- 
insulator transition (see Figs |4j [5] [6|, which coincides 
with the divergence in the "antiferromagnetic" suscep- 
tibility. On the other hand, the DOS for the one-third 
filling (Fig. [7]) does not demonstrate any signatures of 
the Hubbard band formation and the system is very far 



7 




•3 K 
S -6 



-8 
-9 
-[() 



2 electrons, lull 

2 electrons, nn 

3 electrons, lull 
3 electrons, nn 




10 



20 



•!0 



10 



Figure 8. Imaginary part of the self-energy in the Matsubara 
domain plotted as a function of the number of Matsubara 
frequency in the three-orbital model Q with one-half and 
one-third fillings at U = 2.6i, J = U/A, T = 0.05i. The half- 
filled curve shows an insulating diverging behavior, whereas 
the curve for a one-third filling reaches a finite value with a 
proximity to zero. The inset shows the curves for one-third 
filling zoomed at first Matsubara frequencies. 
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Figure 9. Effective mass as a function of temperature in the 
one-third filled three-orbital model (2} at U = 2.6t, J = U/A 
for the density-density and rotationally invariant Coulomb 
interactions. Inset: Imaginary part of the self-energy at uj = 
as a function of temperature for the same parameters. 



from Mott transition. We can attribute this effect to 
a large reduction of the effective Coulomb interaction 
strength from U e ff = U + 2J'm the three electron case 



to U t 
mode 



U — 3 J for two electrons in the three orbital 



However, it does not mean that the one-third filled 
system is not correlated. As the main measure of the 
strength of correlations, one can consider the local self- 
energy function S w . Figure [8] depicts ImE in Matsub- 
ara domain. For the half filling, the behavior is typical 
for Mott materials: X tends to diverge near Fermi level. 
For one-third filling, the system is metallic with a big 
damping due to a strong scattering on spin fluctuations. 
Larger absolute value of the self-energy at the Fermi level, 
obtained from the extrapolation from the first Matsub- 
ara frequencies (inset of Fig. |9| indicates less coherent 
quasiparticles in the "density-density" model while both 



curves with a decrease of temperature approach zero as 
it is e xpect ed for the Fermi liquid and "frozen moment" 
phases™. 

The main curve on Fig. [9] shows that while the in- 
troduction of the full rotationally-invariant interaction 
brings a minor difference to the density of states (on Fig. 
[7]) and the absolute value of the self-energy (in the inset) 
the effective mass is much more susceptible to the par- 
ticular form of the interaction. This is a manifestation 
of the strong correlations in the system. For the density- 
density type of model for J = U/A with the decrease 
in temperature the effective mass is « 2.5 times smaller 
than for the model with the rotationally-invariant inter- 
action. Several mechanisms of this effect are possible. 
First, the stronger enhancement of an effective mass in 
the model with the rotationally invariant interaction can 
be attributed to the 1.5 times larger degree of degener- 
acy of the ground state of the atomic modeF^. Beside 
that, in the described models correlated electrons inter- 
act with effective spin fluctuations. The Ising-like inter- 
action produces fluctuations only for the z-component of 
a spin, whereas the Heisenberg type of interaction im- 
poses quantum fluctuations for all spin components. In 
the simplified description, where the interaction between 
electrons is replaced with the interaction with the fluctu- 
ating classical field, one expects that the electronic self- 
energy behaves as the square of the field fluctuations. 
The latter in general are bigger for the Heisenberg than 
for the Ising type of fluctuations. 

In the standard single-band Hubbard model, s = 1/2 
single-particle local momenta appears because of the U - 
term in the Hamiltonian. The peculiarity of multiorbital 
physics is a formation of a many-electron multiplet pic- 
ture. Taking into account the exchange term J with full 
rotational invariance is of a principal importance here, 
because it is responsible for the Hund's rule validity for 
many-electron configurations. The difference between 
the two cases can be attributed to the mechanism of lo- 
cal moment formation. For the U term, the formation 
of local momenta is accompanied with the localization 
of electrons; Hubbard bands appear due to the charge 
fluctuations. On the other hand, Hund's mechanism of 
a local moment formation does not require localization; 
momenta can appear in a well-developed metallic state. 
One can refer to the t erminology of 'Hund's metal' intro- 
duced previousl y! 48149 1 . 

Clearly, the difference between the rotationally invari- 
ant and Ising-type interaction is more pronounced for 
the Hund's metals. Indeed, charge fluctuations almost 
do not contribute to the electronic Green's function, the 
U term is not of importance, and correlation effects are 
related with spin fluctuations. In this situation, a par- 
ticular form of the exchange part of Hamiltonian must 
affect the result strongly. Larger effect of the symme- 
try of Hamiltonian on the transition temperature for the 
one-third filled system (Fig. [3]) and bigger enhancement 
of effective mass (Fig. |9| confirm this statement. 



VI. CONCLUSION 

We conclude that the magnetism of degenerate orbital 
Hubbard-like models away from the half filling is domi- 
nated by the Hund's rule effects. In this case it is of cru- 
cial importance to keep the proper spin-rotational sym- 
metry of interaction Hamiltonian for the estimation of 
the critical magnetic temperatures and effective masses. 
The corresponding single-particle properties are normally 
metallic while effective Coulomb interaction is small. In 
the half-filled case, this situation is opposite: the domi- 
nating physics is normally related with the large effective 
Coulomb interactions and so the antiferromagnetic tran- 
sition temperature is not very sensitive to the particular 
symmetry of the Hund's rule term. The single particle 
spectrum in this case can be quite sensitive to the spin- 
flip terms due to the Kondo-like effects. 
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Appendix A: Dual fermion perturbation theory 

In this appendix we provide the analytical formalism to 
explain the appearance of the central peak in the density 
of states of the half-filled model with a rotationally in- 
variant interaction as compared to the "density-density" 
counterpart. It was previously noted by Pruschke and 
Bulla^ that the particular structure of the ground state 
of all spins aligned in the same direction in the "density- 
density" model does not allow for the appearance of 
the central peak. This description is clearly valid in 
the Kondo limit, however for the system in the strongly 
correlated regime the argumentation requires some ex- 
tra justification. To proceed with this case, we employ 
the "hybridization-like" expansion in auxiliary fermionic 
variables (dual fermions) around an exact solution of the 



atomic problem^] The vertex functions of the atomic 
problem are used in the diagrammatic expansion. The 
appearance of the peak is connected to the divergence 
in the dual fermion self energy. It will be shown that in 
the model with a rotationally invariant interaction this 
divergence is guaranteed by the structure of the vertex 
function. 

One important circumstance is that for a degener- 
ate atomic ground state such a theory cannot be con- 
structed straightforwardly, since the degeneracy, in gen- 
eral, breaks the perturbation series. In a work by Hafer- 
mann et aiP^ the dual fermion perturbation theory is 
written for a symmetric single-band Anderson model (see 
figure 4 of [50 ) . It has been shown that the Kondo-like 
resonance is correctly reproduced within this approach 
only if the reference system ("isolated atom") is in fact 
not just a single impurity atom, but a cluster includ- 
ing the impurity and at least one bath site ("superper- 
turbation theory"). The superperturbation theory does 
not suffer from the degeneracy issue, because the ground 
state of such a cluster is a singlet. 

An alternative way proposed in [40] is to deal with a 
single isolated atom, and to break the symmetry of its 
ground state with a small external field, so that the per- 
turbation series is constructed around a polarized atomic 
solution (a similar approach was proposed by Logarpif). 
Such a symmetry-broken version of the theory makes a 
qualitative analysis feasible. 

For simplicity, we consider the paramagnetic phase 
only (that is, put h — in (|6|) and write the formulae 
for the zero temperature. The impurity action S lrap is 
split into a sum of an atomic action corresponding to the 
Hamiltonian ^ and a hybridization term (a summation 
over orbital and spin indices a is supposed): 



S imp = 5 at - dt h a n ta 



+ 00 



dt dt' c* ta (A a (t - - h a S(t - t')) tva, (Al) 



where A a (t — t') is the self-consistent hybridization func- 
tion. The integrals are taken over the axes of real times t. 
An infinitesimal field h a is introduced to break the sym- 
metry of the atomic problem; the result will be averaged 
over all possible orientations of h a . 

Following the procedure of the dual transformation^, 
a new set of fermionic (dual) variables f* a , f ea is then in- 
troduced through a Hubbard-Stratonovich decoupling of 
the last term in jAlj ): c*(A - h)c -> cf* + c*f - /*(A - 
h)^ 1 f. An average value, represented by a path integral 
over c*,c with a kernel exp(iS* lmp ) now can be rewrit- 
ten in the introduced mixed representation as an integral 
over both c*,c and new variables /*,/. Moreover, it is 
possible to switch to the pure dual representation by in- 
tegrating out only c*,c variables. This leads to a dual 
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action with, in general, nonlinearities of all even orders: 

+00 

5 duai = -J2 J de (G° a (e))- l f;j ea +V d [r,f] (A2) 



G°(e) = 



1 



(fc(e)- 1 - A a (e) 



-3a(e) 



(A3) 



Functions <? Q (e) are one-particle Green's functions of the 
atomic problem: g a (t — if) = —i(^f\Tct a c^ a \^) , where 
T is the time-ordering operator and g a (e) is the Fourier 
transform of the g a (t — t'), g a (s) is a-diagonal for the 
models under consideration. Coefficients of quartic and 
higher terms of dual potential V d [f*, f] are proportional 
to irreducible vertex parts of the atomic problem (with- 
out fermionic legs). 

It can be shown that a partition function and any ob- 
servable of the physical electrons is rigorously related to 
the dual counterparts. In particular, for the impurity 
Green's function we have 



The dual fermion representation allows to construct 
a conventional diagrammatic perturbation theory for £ 
using G°(e) as bare lines on the diagrams, and the ir- 
reducible vertex parts as frequency-dependent expan- 
sion parameters. All standard statements, including 
the Wick's theorem, the linked-cluster theorem and the 
Dyson equation hold for the dual perturbation theory. 

One can understand the dual perturbation theory for 
an impurity model as an alternative formulation of the 
strong coupling expansion. Indeed, the zeroth order of 
the theory coincides with the Hubbard-I approximation, 
but the dual transformation gives us Wick's theorem and 
therefore the ability to develop a simple diagrammatic 
perturbation theory. The conventional formulation of the 
hybridization expansion is lacking the Wick's theorem, 
because the isolated atom is an essentially nonlinear sys- 
tem. Thus the achievement of Wick's theorem cannot 
come at no price, and the price for it is the emergence of 
the higher nonlinear terms in the dual action. 

We are going to estimate the lowest-order contribution 
to the dual self-energy ("Hartree" diagram). 



G a (eY 



(A4) 



+ OO 



S Q (e) 



2tt 



d^ [ X^;e,e')^) (A5) 



The dual self-energy £ a (e) is defined through a stan- 
dard Dyson equation for the dual particles. Just like in 
the case of physical electrons, the low-frequency behav- 
ior of the self-energy for the dual fermions reflects the 
occurrence of an in sulat ing gap. 

We can see from (A4| that S(0) — > 00 leads to a metal 
DOS obeying the Friedel sum rule, and S a (0) —> leads 
to the Hubbard-I solution with an insulating gap (which 
may be of zero width). This fact is somewhat opposite 
with respect to the analogous statement for the physi- 
cal self-energy. But it comes as no surprise, since the 
physical S tt (e) is an essential building block in the weak 
coupling theory and, as it mentioned above, we are deal- 
ing with the strong coupling expansion. 



Note that this diagram equals zero for the dual series 
built on top of a DMFT solution. However, here we con- 
struct an expansion on top of the atomic problem, and 
(A5) does not vanish. Moreover, we will discuss in what 



case it diverges at the Fermi level. In accordance with the 
discussion above, such a divergence indicates a pinning 
of the metallic DOS. 

Let us analyze expressions for the two-particle Green's 
(4) 

function X1234 an d the corresponding vertex function 
M234- F° r rea l frequencies at T = there exists a spec- 
tral representation for X1234 m analogy to expressions 
given in Refs. 1521 and I5U1 



Xi234( £ i > £ 2; £3, £4) = s S n ( n ) X^*l^ n i l fc ) ( fc l^ n 2 \ l ) ( l \On 3 \m) (m|d n J*) faimiznt , zn 2 , zn 3 , z n J 



/,7 m 



(A6) 



'k-lni 



(z ll z 2 ,z 3 ,z 4 ) = -?; 



z\ + Em — E k 4- iS z\ + z 2 + Ey — Ei + iS z\ + z 2 + z 3 + E^ - E m + id 



(A7) 



r 



where \&) is the ground state, indices k,l,m denote the 
eigenstates of the Hamiltonian ^ and n denotes the 
permutation of the operators, which happens due to the 
time-ordering. The sum over n is thus a sum over all 
24 permutations of operators O = {c\, C2, Cg", C4 }. The 



whole frequencies dependence of the vertex is contained 
in the spectral kernel 4>Mm- Arguments z are defined as 
{ei, £2,-63, -£4}- 

The symmetry breaking of the atomic problem by an 
infinitesimal external field makes the ground state $ 
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unique and so no summation over \I> is performed. The 
k, m states differ from VP in the number of particles, and 
have therefore a different energy (see Tables |lj |n| . Con- 
trary, Ei can coincide or be very close to Eq, (if I coincides 
with or another ground state of the unperturbed atom, 
respectively). This "singular" case is particularly inter- 
esting to us, because a resonant kernel [z\ + Z2+iS)~ 1 
may turn into (e — e' ± i<5) -1 in the formula (A5) and 
result in a logarithmic divergence at e -> as soon as 
Gfl(e') goes to a finite value at the Fermi level. 

The multiplet structure of the atomic ground state is 
of crucial importance here. First of all, if the ground 
state of the atomic problem is unique, a singular term 
in the two particle Green's function will inevitably cor- 
respond to |Z) = \^f). But such singular terms cancel 
with the one-particle reducible part and thus do not en- 
ter the vertex. From this we conclude that singular terms 
in r^2 3 4 appear, in general, for atoms with a degenerate 
ground state. Furthermore, not all of the states of the 
multiplet contribute to the singular part, since \l) is con- 
nected to |^) by a pair of creation-annihilation operators 
On~i * On 2 . It means that a contributing term is produced 
only when total spin projections of states \l) and \^) dif- 
fer by ±1. 

A resonant term (e — e' ± iS)^ 1 in the integrand of 



(A5) results in the logarithmic divergence of S(e = 0) 



and, consequently, in the pinning of DOS at Fermi, when 
G°J0) is finite. The latter is true if both A /3 (0) and 
5/3(0) are finite. When the first part of the condition 
means the bare density of states being conductive, the 
last one is satisfied when the particle-hole symmetry is 
broken in the corresponding /3-mode. But as we consider 
solutions with broken symmetries, there must be at least 
one combination of indices (3 and | which gives gpiO) ^ 
0. 

Now it is easy to see that an impurity model with the 
density-density interaction does not produce a pinned so- 
lution, whereas rotationally-invariant does. The ground 
states of these models are two-fold degenerate (see tables 
I and II). Components of the doublets have the same 
orbital moment projections but their spin moment pro- 
jections differ by AS Z — 1 for the 2-band model and by 
AS Z = 3/2 for the 3-band model. There is no such a pair 
of operators c a cj which could connect the components of 
the doublets, because the differences AS Z are too big. As 
a consequence, singular terms in are not produced 

by density-density interactions in the electron-hole sym- 
metric case and the resulting impurity DOS turns to be 
insulating within the chosen approximation for the dual 
self-energy. 
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